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1. INTRODUCTION 



conditions at the surface are purely diffusive. 



In recent years, the behavior of quantum gases 
has increased interest. In particular, this is related 
to the development of experimental procedure for 
producing and studying quantum gases at extremely 
low temperatures [l|. In the majority of papers, bulk 
properties of quantum gases have been studied [2[ and 
[3fl. At the same time, it is obvious that it is important to 
take boundary effects on the properties of such systems 
into account. In particular, such a phenomenon as the 
temperature jump at the interface between a gas and 
a condensed (in particular, solid) body in the presence 
of a heat flux normal to the surface is important. Such 
a temperature jump is frequently called the Kapitsa 
temperature jump [4j. We also mention a paper where 
the problem of boundary conditions for the motion of 
super liquids was considered |5|. 

The problem of a temperature jump in a quantum 
Fermi gas was studied in [6|, where an analytic solution 
for an arbitrary degree of gas degeneracy was obtained. 
A similar problem for a Bose gas was considered in 
[7], where the gas was assumed to be nondegenerate, 
i.e., it was assumed that there was no Bose — Einstein 
condensate. 

The present paper is devoted to solving the problem of 
a temperature jump in a degenerate Bose gas analytically. 
The presence of a Bose — Einstein condensate |8( leads 
to a considerable modification of both the problem 
statement and its solution method. In this kinetic 
equation with a model collision integral is used to 
describe kinetic processes. We assume that the boundary 



'-V/ = /[/], 



2. KINETIC EQUATION 

To describe the gas behavior, we use a kinetic equation 
with a model collision integral analogous to that used to 
describe a classical gas. In this case, we take into account 
the quantum character of the Bose gas and the presence 
of the Bose — Einstein condensate. 

For a rarefied Bose gas, the evolution of the molecule 
distribution function / can be described by the kinetic 
equation 

df <9£, 

— H 

dt dp 

where £ is the kinetic energy of molecules, p is the 
molecule momentum, and /[/] is the collision integral. 

In the case of the kinetic description of a degenerate 
Bose gas, we must take into account that the properties of 
the Bose — Einstein condensate can change as functions 
of the space and time coordinates. In other words, we 
must consider a two-liquid model (more precisely, a two- 
fluid model, because we consider a gas rather than a 
liquid). We let p c = p c (r,t) and u = u(r,i) denote the 
respective density and velocity of the Bose condensate 
Einstein. We then have the expressions 

„ Pr.i' 2 

J = Pc", 



Tlik = PcUiUk- 

for the densities j and Q of the respective mass and 
energy fluxes and for the momentum flux tensor lT^j. of 



2 



the Bose condensate — Einstein (under the assumption 
that the chemical potential is zero). 

The conservation laws for the number of particles, 
energy, and momentum require that the relations 



dp c 
dt 

dE c 

dt 

d{p c u) 
dt 



vj = - J i[f}dn B , 



vq = - 



vn 



2m 



I\f\d£li 



pi[f]dn B 



be satisfied, where s is the molecule spin 
(2s + l)d 3 p 



PcU 



(27Tfi) 3 ' 

where h is the Planck constant. 

In what follows, we are interested in the case of 
stationary motion with small velocities (compared with 
the thermal velocities). We note that for the Bose 
condensate, the quantities Q and , are depend 
nonlinearly on the velocity (they are proportional to 
the respective third and second powers of the velocity) . 
Therefore, in the approximation linear in the velocity 
u , the energy and momentum conservation laws can be 
written as 



P 
2m 



and 



J P i[f]dn B = o. 



According to the Bogolyubov theory, the relation for 



the excitation energy e(p) 



e{p) 



2 2 
U p 



P 

2m 



1/2 



(1) 



where 



4nh 2 an 



1/2 



Here a is the scattering length for gas molecules, 
n is the concentration, m is the mass, and p is the 
momentum of the gas molecule, holds for a weakly 
interacting Bose gas. The parameter a characterizes the 
interaction force of gas molecules and can be assumed to 
be small for a weakly interacting gas. 



The relation u 2 <C kT/m , where k is the Boltzmann 
constant and T is the gas temperature, holds for 
sufficiently small a . In this case, we can neglect the first 
term in the brackets in (1). The expression for the energy 
£(p) takes the same form as in the case of noninteracting 
molecules: 



£(p) 



p 

2m 



We now consider the widely used kinetic equation 
in the Boltzmann — Krook — Welander form with 
the molecule collision rate proportional to the molecule 
velocity 



dj_ 
dt 



+ (W)/ = v w(f* M - /). 



(2) 



Here, / is the distribution function, v is the molecule 
velocity, w = |v — vo| , vo is the mean-mass gas velocity, 
f^j is the Maxwell distribution function, 



A I 



( 771 



3/2 



exp 



2fcT» 



(v- 



and v is a model parameter corresponding to the inverse 
mean free path / of a molecule, vq ~ l/l . 

The parameters in the formula for f^j , namely, ti* , 
and u* , can be determined from the conservation 
laws for the number of molecules, momentum, and energy 



wf dfl 



wfijdnM, 



wvf dQ M = / wwf* M d£l 



(3a) 



(36) 



w— (v - u) 2 fdn M 



where 



dn 



w-{v-xi) 2 f* M dn M , (3c) 



d A v. 



The conservation law for the number of particles in 
the normal state is inapplicable because the transition of 
particles to the Bose — Einstein condensate can occur. 
As mentioned above, the effect of the condensate on 
the energy and momentum conservation laws can be 
neglected in the approximation linear in u . 

We note that Eq. (2) corresponds to the assumption 
that the mean free path of molecules is constant (it is 
independent of molecule velocities) . It hence follows that 
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Eq. (2) to a greater extent corresponds to the model in 
which the molecules are regarded as solid spheres. 

We consider a generalization of Eq. (2) to the case 
of a degenerate Bose gas. We assume that the general 
structure of Eq. (2) is preserved, but by a function , 
we must mean the Bose — Einstein distribution (Bosean) 
with a zero chemical potential 



f*B = 



( m I 



u,)" - 1 



Here, the parameters T* and u* , are determined by 
the second and third conditions in (3). In this case, we 
have 



2s + l_^ 3 
(2tt/i) 3 



instead dfl]\j . 

We assume that the mass velocity of the gas is much 
less than the mean thermal velocity of the molecules 
and the typical temperature variations along the mean 
free path of molecules are small compared with the gas 
temperature. The problem can be linearized under these 
assumptions. 

We seek the distribution function in the form 



f = f s B (v)+v(t,r,v)g(v), 



where 



i 



exp((3 s v 2 ) - 1 



(p is a new unknown function, T s is the surface 
temperature, and 



We introduce the notation 



(v-u*) 5 



2kT, 



Taking this notation into account, we have 



/1(C) 



exp(e*) — 1 



exp(C 2 ) - 1 ' 



9(C) 



exp(C 2 ) 



(exp(C 2 ) - If 



We linearize the local Bose — Einstein distribution fg , 
passing to dimensionless quantities. We note that 

£ » = ^[/3,(v-u,) 2 ]=^[(C-W) 2 ], 



T, 



where 



Taking = T s + 5T S into account, we obtain 
e» = C 2 - 5 -^C 2 - 2CW, 



whence we find 



6e„ 



o 5T* 

-2CW-C 2 — 



where 

Consequently, 

JB = JB 

or 

r B = f B +9(c) 



dfh 



5e*, 



2CW + C' 



We note that the quantity w in Eq. (2) can be replaced 
with v in the approximation under consideration. We 
introduce the dimensionless quantities t* = tVo/j3 s and 
r* = ri>Q and omit the asterisks on these quantities 
below. It is now clear that Eq. (2) (in the dimensionless 
variables) becomes 



g + (CV) V = 2C(CW) + C^ 



Op. (4) 



The parameters of this equation can be found from the 
momentum and energy conservation laws (relations (3)), 
which now become 



{L(f)C 2 g(C) d 3 C = 0, 



where 



j (L<p)Cg(C) d 3 C = 0, 



Crji 

Lip = 2CCW + C 3 — - - C<p. 

J- « 
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From this system, we find 
3 



W = 



6T t 1 



T s 4irg 2 



V3C.g(C)C(i d C, 



J vg(C)C 3 d 3 C, 



where 



oo 

9n = J C 5+n g(C)dC, 



n = 0,1,2, 



In this case, we have 



oo oc 

5o = J g(C)C 5 dC = J 



cxp(C 2 )C 5 dC 

{ (cxp(C2)-l) 2 



-2 J Cln(l -cxp(-C 2 ))dC 



1.64493, 



temperatures. The quantity A is called the temperature 
jump (the Kapitsa temperature jump in the case of low 
temperatures). 

The problem is to find AT as a function of the heat 
flux Q . Taking into account that the problem is linear, 
we can write 

AT = T Q Q. 

The quantity Tq is called the coefficient of the 
temperature jump. Another notation can be used: 

Q = i?AT, 

where R is called the Kapitsa resistance. 

Taking into account that the problem is stationary and 
that the function ip is independent of the coordinates 
and z , we simplify Eq. (5) : 



3i = J g(C)C 6 dC = 2.22912, 
o 

oo 

5 2 = J g(C)C 7 dC = 3.60617. 

o 

We represent Eq. (4) in the standard form 

dip 



dt 



+ (CVV + C^,r,C) 



= £ J k(C, C>(i, r, C) C'3(C')d 3 C", (5) 



rfle 



fc(C,C) = — CC' + — c 2 c' 2 . 
9a 32 



3. PROBLEM STATEMENT 

A degenerate Bose gas occupies the half-space x > 
above the plane surface in the problem under 
consideration. A heat flux Q normal to the surface 
is maintained in the gas. We let T denote the gas 
temperature far from the surface. The quantity To differs 
from the surface temperature T s if there is a heat flux. 
We let A = T — T s be the difference between these 



± oo 

\ 1 1 k(n,C;n',C')<p(x,»',C')<Kl(C), (6) 



where 



-1 o 



/i = — , dn(C') = g{C')C lZ dfi'dc', 



k(c, n\ a, m') = —c^Cfi' + —c 2 c 2 . 



go 



52 



It is easy to verify that Eq. (6) has the particular 
solutions 

ipi = nC, tp 2 = C 2 . 

Consequently, the function 

ip as {x,n,C) = ByiC + e T C 2 , (7) 

where is proportional to the heat flux, is an asymptotic 
distribution function (as x — > +oo ). 

Assuming that the reflection of the molecules from the 
wall is purely diffusive, we now formulate the boundary 
conditions 



ip(0,fj,,C)=0, 0</i<l, 



<p(x,n,C) 



(8) 



= ip a3 (x,fj,,C) +o(l), x -> +oo, -\<fi< 0. (9) 
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4. REDUCTION TO THE ONE VELOCITY 
PROBLEM AND THE SEPARATION OF 
VARIABLES 

We seek the solution of problem (6)-(9) in the form 

<p(x, fi, C) = Ch x {x, fi) + C 2 h 2 (x, fi). (10) 
We obtain the system of equations 
\i— + hi.{x,ii) = 



1 1 

(i J n'hi{x,n')dn' + ^-/j, J n'h 2 {x, (i')dn', 



dh2 , u ( \ 



9i 
252 



J hi(x, fj,')dfi' 



+ / h 2 (x, fi )d[i . 



We represent this system of equations with respect to 
the column vector 



h(x,fj,) 
in the vector form 



hi(x,n) 
h,2{x,n) 



i 

M <9z + h ( x,fJ ') = \ J K(^,fi')h(x,fi')dfj,', (11) 

-l 



where K{n,^!) is the kernel of Eq. (11), 



3/U/li' 3 — HfJ,' ^ 



9i 

92 



9o 
1 



Using (10), we transform boundary conditions (8) and 
(9) into 



h(0,n) = 0, 0<fi<l, = 



(12) 



h(x,n) = 

h as (x,n)+o(l), x ->■ +oo, -1 < fi < 0, (13) 



where h as (x,/j l ) is the asymptotic distribution function 



We note that the equation kernel can be represented 
in the form 

K(n,n') = Ko + Snn'Ki, 



K 



The general method for the separation of variables 
yields the expression 






o" 




1 


9i ' 


9i 


1 






5o 


. 92 











hr,{x,n) = exp(--)$(r?,/i). 



(14) 



Substituting (14) in (11), we obtain the characteristic 
equation 



{ri-n)$(ri,n) = |^ n (0) (??) + -^i^ (1) (»7) J 



where 



1 

(r?) = J ^ k <P(i]^)dfi, k = 0,l. 



It is obvious, that 



( V ) = r 1 (E-K )n^(v), 



where E is the unit matrix of the second order. 
We obtain now the characteristic equation 



(15) 



where 



n{rj) = n (0) (rj) = J $(r), fi) dfi, 



D( m ) = K + SftriK^E - K Q ) = 



3g/irj 

92 



5 = 1- 



5o 92 



= 0.96288. 
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For rj e (—1,1) , we use Eqs. (15) in the class of 
generalized functions [ill to find the eigenvectors 

®{V,V) = F{ri,^)n{ri) 
of the continuous spectrum, where 

F(V, M) = \riD(jii))P^— + A(v)S(v - n), 

2 n-n 



n{rj) 



(16) 



F(rj,/j) is the eigen matrix-function, the symbol Px^ 1 
denotes the principal value of the integral in the 
integration of x~ x , 6(x) is the delta function, and A.(z) 
is the dispersion matrix, 

l 

, , . „ z f D(uz) , 
A(z) = E+- / -^dfi. 

2 J fJ. - z 

-l 

The dispersion matrix has following elements 
A n (z) = l + 3gz 2 X (z) = Ai(z), 

A 12 (z) =0, 



A 21 (z) = ^Ao(z) - ^, 

52 52 



A 22 (z) = A (z). 
We represent the dispersion matrix in the explicit form 
A{z) = X (z)D(z 2 ) + D , 

where 



D 



1 

5i 







52 



A (z) = 1 + 



du 

2 J u — z 
l 



Its determinant, called the dispersion function, is given 



5. EIGENVECTORS OF THE DISCRETE 
SPECTRUM 



By definition (7|, the set of zeros of the dispersion 
function is called the discrete spectrum. 

The function Xq{z) has a single double zero at the 
point at infinity Z{ — oo . Two eigensolutions of Eq. 
(11) correspond to this zero (they coincide with the 
eigenvectors of the characteristic equation): 



ho(x,fi) = /i 



hx(x,fi) 



We find zeros of Xi(z) outside the cut [—1,1] . It is 
obvious that Ai(oo) = 1 - 5 = 0.03712 > . We take the 
contour j e (see Fig. 1) that encompasses the cut [— 1, 1] 
at a distance e from it such that there are no zeros 
of the function Xi(z) inside the contour. The number 
N of zeros of Xi(z) outside the contour is equal to the 
increment of the argument of this function according to 
the argument principle 12 1, i.e., 



7V=(2 7 r)- 1 A 7e argA 1 (z), 

where the symbol A 7e / means the increment of / along 
7 e . Passing to the limit as e — > in this equality, we 
obtain 

iV^A^arg^ 



2tt 



where Xf(fi) are the boundary values of X±(z) on the 
interval (—1,1) from above and from below, 

XfifJ.) = Ai(/i) ± Siirgn 3 . 

Taking into account that Reu + (fi) = w(/i) is an even 
function and Imw + (/Lt) is odd, we have 



N = -A (M axgG(ji), 

7T 



KM' 



Let 9i(p) = argAj (/x) be the principle value of the 
argument specified by the condition Q\ (0) = . Because 



by 



A+( M ) = Ar(/i), |a+(m)| = |ArO0l, 

we have arg G(/i) = 25i(/i) and therefore 



X(z) = det A(z) = X {z)Xi(z). 
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has a nonzero solution because 

detA(r7o) = A(?7o) = 0. 

We represent Eq. (17) in the form of two scalar 
equations 



Ai(»7o)rai(»?o) = °> 



(18) 




Fig. 1. 



Fig. 1. Dependence of the temperature jump coefficient on 
the parameter 7: the specular reflection coefficient is q = 0.3 
for curve 1, q — 0.5 for curve 2, and q = 0.8 for curve 3. 



It is easy to see that the angle 0i(/z) on the cut 
[0, 1] has an increment equal to tt ; consequently, we 
have N — 2 (the number of zeros is equal to two). We 
let ±7/o (r/o = 1.27573) denote these zeros. In view of 
the equality Xi(z) = Xi(z) , these zeros are real. Two 
eigensolutions h± Vo (x, /j) , where 



h no {x,n) = $(?7o,/x)exp( ), 



Vo 



2 r/o - n 

correspond to them. 

We note that the homogeneous equation 



(17a) 



(176) 



9i 

fj2 



[AoM - l]ni(r? ) + A (?7o)n 2 (?7 ) = 0. (19) 



In view of the condition Ai(?7o) = , it follows from 
Eq. (18) that the upper element ni(r]o) is arbitrary and 
nonzero. From Eq. (19), we now hnd the lower element 
of the vector n{rj ) , 

n 2 {vo) = ~— (l - x 7 , )ni(m)- 
52 V A (?7 )y 

From the equation 

Ai(r/ ) = 1 + 3.gr/oAo(r7o) = 

we find 



1 



Hence, the vector n(r/o) has the form 

1 



n(r? ) 



9i 
92 



(l + 3 OTo 2 ) 



n\(r] ). 



We substitute this vector in the second condition in 
(17) and obtain 



for 



ni(??o) 



1 



9i, 

92 



-2A (r/o). 



6. DECOMPOSITION OF THE SOLUTION INTO 
EIGENVECTORS: THE RIEMANN HILBERT 
BOUNDARY VALUE PROBLEM 

We show that the solution of boundary value 
problem (11)— (13) can be represented in the form of a 
decomposition, namely, 



h(x, n) = h as {x, (j,) + A a exp( )$(r/o, A»)+ 

Vo 



exp(— )F{r 1 , l i)A(r])dr], (20) 
V 



where the unknowns are the constants Et and Aq and 
a vector function A(jf) with the elements Aj(j)),j = 
1,2,3. Decomposition (20) can be represented in the 
form 



h(x, n) = h as (x, + A Q exp( )*(??o, £*)+ 



-exp( — )A(p)A(jj)0 + (ii)- 
A* 



-J / eM--)vD( m )A(v) — , (21) 
2 7 V rj-fj, 

o 



where 



_( At ) = l lM e(0,l); e + {n) = 0,ni(-l,0). 



Let's lead reduction to the diagonal form of the 
problem (25). For this purpose we will lead to a diagonal 
kind the matrix 



P{z) 



3gz 2 
ffi 



o 

Ao(z) 



The matrix bringing the matrix -P(z) to a diagonal 
kind, is that 



S = 



-1 

1 a 
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It is obvious that 



detS = l, S- 1 = 



91 

92 



1 



-1 



S- 1 P(z)S = fi(*) = diag{^, X (z)}. 

We first solve the homogeneous boundary value 
problem 

P+(ji)X+(n)=p-(ii)X-(it), < a* < 1, (25) 



Substituting decomposition (21) in boundary condition 
(12), we obtain a singular integral equation with the 
Cauchy kernel [12j 

h as (0,n) + A $(?7o,At) + A(n)A( f i) + 



~ [r l D( m )A(r,)-^- = 0, < /x < 1. (22) 
2 7 ?7 - m 

o 



We introduce an auxiliary vector function 
l 



N(z) = \ f v D(zr,)A(r,)-^- 
2 J r)-z 
o 



corresponding to (24). 

Clearly, that it is necessary to search for matrix X(z) 
in the form 

X(z) = 5diag{c/!(z), U 2 (z)}s- 1 . 

The method for solving such problems was developed 
in [jj; therefore, we give the solution of problem (25) 
without a derivation, 



X(z) = 



(23) 

where 



9x 
.92 ' 



Ui(z) 
[U {z)-U 1 {z)\ Uo(z) 



and the matrix 



U a (z) = zexp(-V (z)), 



P{z)^k{z)D- 1 (z 2 ). 



U 1 (z) = zeM-Vi(z)), 



Using the boundary values N(z) , A(z) and P(z) 
and the corresponding Sokhotski formulas, we reduce 
Eq. (22) to an inhomogeneous vector Riemann — Hilbert 
boundary value problem 



V Q (z) = 



1 f Co (u) du 



ir j u — z 
o 



P+(/i)[iV+(/i) + h as (0,(J,) + A $(r; , M)] = 

p-(jx)[N-{fi) + h as (0,fi)+Ao^(rio,iJ,)], < a* < 1. 

(24) 



Vi(z) 



i f Ci ( M ) ^ 



7T J M — 2 




Co(«) = e («)-7r, Co(0) = -tt, Co(l-0) = 0, 
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Cl(«)=0l(u)-7T, Cl(0) = -7T, Cl(l-0)=0. 

In this case, the angles have the form 

2\ {u) 



7T 



Co(u) = - - - arctg- 

Z TTU 



a / \ n 2Ai(u) 

ClW = -77 - arCt S "q a"" 



We now return to the solution of inhomogeneous 
problem (24). Substituting the matrix P + {^) found 
using (25) in (24), we obtain the problem of determining 
an analytic vector function from its jump 



X + {ii) Ar+(/x) + /i a5 (0, M ) + Ao$(77o,M) = 



= X-(ji)\ [N-(ii) + h as {0,ri + A $(ri ,iJ,)\. (26) 
We note that 



X~\z) 



1 



9i 



Ui(z) 
U U 

\Un(z) UUz)J 



L 92\U {z) U^z)} Uo(z) J 
N(z) 



Ni(z) 
N 2 (z) 



The asymptotic behavior of these functions in a 
neighborhood of the point at infinity can be described 
as 



X~\z) 



N(z) = 



1 



z- 2 z- 1 







(0) 



+ o(l), 



z — > 00, 



00, 



where 



< 5 = -| Jv 2 Mv)dv 


according to (23). 

Taking the behavior of X~ x (z) and N(z) at finite 
points of the complex plane and in the neighborhood of 
the point at infinity into account, we obtain a general 
solution of problem (26) 

N(z) = -h as (0, z) - A $(t7o, z)+ 



+X(z) 
or, in scalar form, 
N^z) = -Bz- 





■c x - 


1 


■di- 






- - 


+ 

Vo - z 


d 2 - 





A z 



+ U 1 {z) 



Ci 



di 



Vo- z 



(27) 



at f \ ,9i ^oVo . 
N 2 (z) = -e T H h 



92 Vo - z 



9i 

92 



U {z)-U 1 {z) 



di 



Vo- z 



d2Uo{z) 
Vo- z 



• (28) 



Here C\ , d\ and d 2 are constants that can be found 
from the solvability conditions for the boundary value 
problem. 

We note that 

U (z) = z - Vo" 1 +o(l), z -> 00, 
J7i(z) = z - Vf 1 +o(l), z ->■ 00, 



where 



1 

V/ _1) = ~ J Ci(w) du = 0.84188, 



V 



(-1) 



(u) du = 0.71045, 



Eliminating the pole of the function N\(z) given by 
equality (27) at the point at infinity, we obtain C\ = B . 
Equating the limits on the right and on the left at the 
point z = 00 in equality (27) for N\(z) and using (23), 
we obtain the equation 



2 



J v 2 Mv) dv = -A + di + BV} 1} . 



(29) 



From the equation N 2 (oo) = 0, where the function 
N 2 (z) is given by equality (28), we find 



e T = -d 2 + B 



9i 
92 



(30) 



Eliminating the poles of the functions N\(z) and 
N 2 (z) at the point 770 , we obtain 



di 



AqVq 
Uiivo) 



10 



and 



, 9i , 9i A oVo 
d 2 = di = -- 



92 92 f/i(lo) 

Taking these relations into account, we can rewrite 
equality (30) as 

AqTJq 



9i 

92 



Ui(vo) 



+ B(V{- 1} - V ( - 1} ) ■ (31) 



From definition (23) of the function N(z) , we have 



2 J rj - z 
o 

whence 



N+(fi)-N^( t i) = 3gni l j, 3 A 1 (fx) 

according to the Sokhotsli formulas. 

Substituting solution (27) in (32), we find 

3 9 . , 2 



(32) 



1 

'2id 



V Tin J 



rial rj Voiv-Vo) 



Uf(ri)-Ur(ri) ■ (33) 



Integrating equality (33) over 77 from to 1 , we 
obtain 



2 

where 



TfA 1 {r,)dT l = (d + -) J(0) - -J(ryo), (34) 
V t/q J Tin 



'/() 



J(z) = J-/^-^dr. 
w 2vrz / t-Z 



We give the integral representation of the function 
Ui(z) without a derivation, 



1 [U+(t)-U^(t) 



2iri 



dr, z £[0,1]. 



According to this representation, we have 
and 

J(0) = C/ 1 (0) + ^ 1) . 



Substituting these equalities in (34), we obtain the 
equation 
1 

35 lv 2 A 1 (r i )dr l =(c 1 + d] 

no- 



rfA 1 (ri)dr,= (c 1 + ^)(u 1 (o) + v}- 1) } 



di 
Vo 



(-1) 



(35) 



From Eqs. (29) and (35), we now obtain an equation 
from which we find 

= -bum. 

Substituting the found value Ao in (31), we find that 
the temperature jump is given by 

,9i 



e T = B- 



92 



■n 4- T/f" 1 ) 
T]0 + V l - V 



(36) 



7. TEMPERATURE JUMP 



We express the temperature jump in terms of the heat 
flux, which is only transferred by the normal component 
and is proportional to its velocity 0]. We note that the 
mean velocity of the gas (the normal component and 
Bose condensate) is equal to zero in the gas volume. We 
calculate the heat flux using the formula 

f mv 2 (2s + l)d 3 p 
Q = / /v— 77T-^ — > p = TOV. (37) 



2 (2ttS) 3 

Passing to dimensionless quantities in the integral in 
(37), we obtain 



Q 



(2s + l)m 4 



2(2vrft) 3 ^ 6 



J [f s B + ip{x,C)g{C)]CC 2 d 3 C, 



(2s + \)mk 3 T 3 

Q " 2^h 3 x 



Chi(x, (j,) + C 2 h 2 (x, n) CC 2 g{C)d 3 C. 
Taking into account that 

C = (C/i, C sin 9 cos x, C sin 9 sin x) , 



we obtain 



d 3 C = C 2 dfidCd X , 



Qx{x) 



11 



Qo 



1 oo 



-1 



Ch!(x,n) + C 2 h 2 {x,^) ^C 5 g(C)d^dC, 



where 



Qo = 



(2s + l)mfc 3 T 3 
^3 ■ 



We calculate the inner integral in (38), 



Qx = Qo(gi J (J,h 1 (x,n)dfi + g 2 J \ih2(x,n)dfxj. 

-l -l 

Taking into account that the heat flux is conserved 
(i.e.. Q x — const ), we find 



hence 



B = 



3Qx 



(39) 



2<9o.9i 

Taking (39) into account, we represent (36) in the form 
AT = RQ X , (40) 
where R is the Kapitsa resistance, 



R 



3 r 

2<72 



ir 2 h 3 



(2s + l)mfc 3 T s 2 ' 



R = 0.58514 



= 5.77510 



ir 2 h 3 



(2s + l)™fc 3 T2 



ft 3 



(2s + l)mfc 3 T s 2 ' 



(41) 



Formula (40) is the sought temperature jump (the 
Kapitsa jump) in the degenerate Bose gas. The coefficient 
R of the temperature jump, called the Kapitsa 
resistance, is given by (41). It follows from relation 
(1) that obtained formula (41) is applicable under the 
condition 



decreases. At present, there are no experimental data 
on the Kapitsa resistance for a degenerate Bose gas. 
There exist only data on the Kapitsa resistance for liquid 
helium 4 He 0| , . According to these data, there exists 
a divergence of the Kapitsa resistance as T s — ► in the 
case of liquid helium. The experimental data thus agree 
qualitatively with formula (41). 
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8. ANALYSIS OF RESULTS 



It can be seen from (41) that the Kapitsa resistance 
increases without bound as the temperature T s 



